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(54) Method and apparatus to determine the speed of mobile communications apparatus 



(57) A method tor determining the speed of a mobile 
user, comprising deriving first signals related to the 
phase variations of a number of taps of an equivalent 
channel input response at the input of a channel equal- 
izer in a base station receiver which receives the signal, 

deriving second signals related to the carrier fre- 



quency offset of the signal received in the base sta- 
tion receiver; 

subtracting the second signals from the first signals; 
and deriving from the difference a quantity related 
to the Doppler frequency shift of the received sig- 
nals, and therefore to the speed of the mobile. 
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Description 

Field of the Invention 

s [0001] This invention relates to a method and apparatus for determining the speed of travel of a mobile telecommu- 
nications apparatus, such as in individual mobile telephone, for example in a GSM (Global System for Mobile Com- 
munications) system. 

[0002] A typical wireless mobile telecommunications system 10 is shown in Figure 1A. The system comprises a 
number of users of mobile transceivers 11, 12, 13, 14, a number of base transceiver stations 15, 16, 17, 18, connected 
10 through Base Station Controllers (BSC) 19, 20, 21 to a Gateway Mobile Services Switching Centre (GMSC) 22. The 
GMSC is connected to the appropriate part of Public Switched Telephone Network 23. 

[0003] A mobile telecommunications system is arranged to have spatial cells through which a mobile telephone 
travels. A recent development is to provide a two tier system having both small-area microcells, which provides in- 
creased network capacity by increasing frequency re-use, and large-area macrocells which reduce the number of hand- 
is . offs as a mobile telephone moves from one cell into the next. 

[0004] Each hand-off requires network resources to redirect the call from one base station to the next, and as the 
time between hand-offs decreases the quality of service may be reduced. 

[0005] It would be highly advantageous if fast-moving mobiles could be allocated to macrocells (to reduce the need 
for hand-offs) and slow-mobiles could be allocated to microcells, so that they can operate at a reduced power and 
20 cause significantly less interference, and minimise the number of handoffs. 

Brief Description of the Prior Art 

[0006] Several methods of measuring the speed of a mobile have been suggested, but all have disadvantages. One 
2S common method is based on measurement of the cell sojourn time of the mobile, but errors may be introduced by the 
different hand-off strategies. Another method measures the signal strength at neighbouring base stations, but accurate 
measurements are not easy to make. 

[0007] Three methods utilise an estimation of maximum Doppler frequencies; the first infers the Doppler frequency 
from the level crossing rate, and the second infers the Doppler frequency from the switching rate of a diversity receiver, 
30 but neither can be used for GSM with frequency hopping. A third methods infers the Doppler frequency from the squared 
deviations of the signal envelopes, but this is not always effective in the presence of multipath fading. 
[0008] A GSM receiver provides estimates of the channel impulse response and the carrier offset so that channel 
equalisation can be performed, but the estimated values have not been used for other purposes than estimating the 
transmitted data. 



35 



Summary of the Invention 



[0009] According to the invention a method of determining the speed of a mobile user in a mobile telecommunications 
system comprises the steps of: (a) deriving a first signal related to a phase variation, in a time interval which is a burst 
40 or fraction of a burst, of one or more taps of an equivalent channel impulse response at the input ol a channel equalizer 
in a base station receiver which receives the signal, (b) deriving a second signal related to the carrier frequency offset 
of the signal received in the base station receiver, (c) subtracting the second signal from the first signal, and (d) from 
the difference, deriving a third signal related to the Doppler frequency shift. 

[0010] Also, according to the invention a mobile telecommunications system comprising at least one mobile trans- 
45 mitter/receiver and a base station having a transmitter/receiver, the base station receiver having a channel equalizer, 
there being further provided first means to derive a first signal related to a phase variation, in a time interval which is 
a burst or a fraction of a burst, of one or more taps of an equivalent channel impulse response at the input of the 
equalizer; second means to derive a second signal related to the carrier frequency offset of the signal received in the 
base station receiver; subtraction means to subtract the second signal from the first signal, and calculation means to 
^0 calculate from the difference a third signal related to Doppler frequency shift of a signal received from a mobile trans- 
mitter/receiver. 

[001 1] Figure 1 B illustrates the Doppler effect in a mobile system. A transmitter T x is shown transmitting at a distance 
€ from the mobile, which moves a distance d from A to B in the direction of the arrow M. 

[0012] Even when the mobile is static, there is multi-path propagation from the transmitter Tx to the mobile, so that 
55 several versions of the transmitted signal arrive sequentially at the mobile. The effect of the different time delays is 
that different phase shifts are introduced between the signal wave components, and there is constructive or destructive 
addition, depending on the relative phases. 

[0013] When either the transmitter or the mobile is in motion, then ■dynamic'' multipath propagation occurs in which 
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a continuous change in the length ol every propagation path appears. This causes changes in the relative phase shifts 
between the paths as a function of spatial location. The time variations in the propagation path lengths can be related 
directly to the relative motion of the transmitter and receiver and directly in the Doppler effects that arise. Figure 1 B 
shows a propagation path from the transmitter arriving first at the mobile at point A along a path of length I, and then 
5 arriving at point B. The time differences between the two points isAt. The mobile moved a distance of d between the 
two points where d = vAt. There is an incremental change A/ in the length of the path where A/ = d cos a and a is the 
angle of arrival of the incoming signal with respect to motion of the mobile. The phase change caused is: 

A0= — r— A/ = — » cos a M 

io XX 

[0014] The rate of change of phase due to motion appears as a Doppler frequency shift in each propagation path 
as follows: 

15 1 A0 v r 

[0015] The change in the path length depends on the spatial angle between the wave and the direction of moving. 
Waves arriving from ahead of the mobile have positive Doppler shift, while the ones arriving from behind the mobile 
20 will have a negative Doppler shift. It is clear from equation C that the higher the speed of the mobile the higher the 
Doppler frequency shift, and the maximum Doppler frequency shift for a certain speed appears when a equals zero or 
1 80 degrees. 

f =^ D 

25 ™ X 

[0016] Doppler frequency correlation with the speed of the mobile can be used to estimate the speed of the mobile. 
However, the fundamental problem in a digital cellular system is that, due to the multipath propagation, each channel 
tap is the sum of many rays. Moreover, due to the memory of the modulated signal, each tap of the equivalent channel 
30 at the receiver is a linear combination of more channel taps. 

Brief Description of the Drawings 

[0017] 

35 

Figures 2 to 15 illustrate various calculated values based on the theory of the invention; 
Figure 16 shows the channel equaliser of a base station receiver; and 
Figure 17 gives details of the implementation of the invention. 

Description of Theory of the Invention 

[0018] The performance of a proposed algorithm has been be evaluated for the case of a GSM receiver. The channel 
model assumed both in the analysis and in the simulation is based on the classic assumption ol omni-directional 
antenna and isotropic scattering. The latter assumption is generally appropriate for a base station operating in a mi- 
45 crocell environment, where the base station antenna is usually mounted at the same height as surrounding objects. 

CHANNEL MODEL 

[0019] Consider the transmission of a GMSK signal s(t). Denoting by l\ the /c-th transmitted bit, and by p(t) the 
50 impulse response of the transmit filter or pulse shape, the baseband-equivalent transmitted signal may be expressed 
by the linear model 



ss 
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[0020] In the mobile radio environment, the signal (1) is modified by the time-varying propagation channel, in which 
short-term signal fading, due to multipath propagation, causes frequency selective distortion as explained above. Time 
variations caused by the Doppler effect may be significant even within the duration of one signal burst. 

5 A. Multipath Channel 

[0021] The received signal is the sum of many Doppler shifted, scaled and delayed versions of the transmitted signal, 
and can be expressed as 

10 

y(t)=x(t)+ w(t) (2 ) 

where 



a(/) = |^-r)c(/,r) C fr > (3) 



20 

c(f,x)denotes the time-varying channel impulse response, and w(\) represents additive white Gaussian noise. At a given 
time t the channel can be regarded as a densely tapped delay line with delay index x 



25 Ay 

= X C n (t)S(T-T n ). (4) 
n-0 



30 
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[0022] Under the assumption that the channel can be modelled as clusters of many independent scatterers, for the 
central limit theorem we can assume Gaussian statistics. Therefore, the channel can be characterised by the mean 
and correlation of c(t, x). When, for a given delay x n , the complex gains due to different scatterers have similar amplitudes 
(diffuse scattering), c(t.x) has zero mean and its envelope follows a Rayleigh distribution. In this case, the tap-gain 
process c n (f) are expressed as 



K-l 



c n {t) = P n (t) = £ (5) 



which represents a cluster of incoherent rays arriving with (approximately) equal delay. In (5), the gain factors ^ the 
Doppler shifts f Dk = f D cos <x k , and the phase shifts 9 k may be assumed fixed during very sort time intervals. Differently, 
in the presence of a dominant path (line-of-sight or specular reflection), the cluster can be modelled essentially as a 
single coherent path 



c n (0= Pn e^W+V. 



(6) 



[0023] In this case eft, x) has a non-zero mean and the envelope has a Ric/an distribution. As a first order statistical 
description of the fading, we can decompose c(t,x) into a specular and a diffuse component. The specular component, 
corresponding to the dominant path, is defined by c s (t,x) = E{c(t,x)} and is known as the channel mean. The diffuse 
component is given by c d (t,x) = c(t,x) - E{c(t, x)} and is Rayleigh distributed. 

[0024] Under the assumption of Gaussianity, a sufficient second order statistical description of the process c(t,x) is 
given by the tap gain cross-correlation function R cc (1 1 ,t, ,t 2 ,x 2 ) = Efcfo ,x, ) c*(t 2 ,x 2 )}. In the case of wide sense stationary 
55 uncorrected scatters (WSSUS) we have 



4 



BNSDOCID: <EP 1052B20A1_1_> 



75 



EP 1 052 820 A1 

It is often possible to write the tap gain correlation function in product form 

Rce& T > = R c<0-*cM 



(7) 



(8) 



where 4> c (t) is the delay power density profile and R c (t) is the normalised autocorrelation junction, R c (0) = 1 The 
10 Fourier transform S c (f) of R c (t) is the Doppler spectrum of the channel fading process. Under the assumption of isotropic 
scattering R c (t) is commonly modelled as 



R c (t) = J 0 {2nf D t) 



where J c (x) is the zero-order Bessel function. The resulting fading power spectrum is 



(9) 
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i/ji-{f// D y 
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M 



(10) 
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B. The GSM Channel Model 

[0025] GSM specifies the propagation models related to different propagation environments, to assess the error 
performance of the receiver. For each propagation condition, the GSM wideband multipath channel is described by 
the time-varying impulse response, with a defined set of tap gains c n (t) = p n (r) ei* n0> of a given mean power and delay 
T n . The real and imaginary part of the complex taps c n {t) are modelled as narrowband stationary independent coloured 
Gaussian noise processes, with variance C n 2 / 2. In the presence of a line-of-sight component Re[c„(t)] and lm[c n (t)] 
have a Rician-type Doppler spectrum, and the envelope p n (t) is Rician distributed. In the absence of a dominant path, 
Re[c n (t)] and lm[c n (t)] are zero mean processes with the classical Doppler spectrum, and the envelope follows the 
Rayleigh distribution. 

[0026] In the GSM channels, the Rician Doppler spectrum is specified only for the first tap of the model and rural 
area conditions (RA). In all the other cases, including the model for typical urban area conditions (TU), the stochastic 
process c„(Q are characterised by the classical Doppler spectrum. 
[0027] Using the definition of R cc one has 

N-l N-l 

JUL',.*****,) - 2 Z 3^)0V,)>*(r,-',)*(*-,-r,)- 

1-0 JmO 

[0028] Now, from the independence of the tap gain processes and under the assumption of Rayleigh distribution (E 
{*„(/)} = 0) 



(11) 
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and therefore 
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4r('i.'*.r) = Z* c M<(**)Wt-t.) = 

= *Jt, - * 2 , r) = E e (t, ~ t 3 ) - <D f (r) (12) 

where, according to (9), 

10 

"o<V'*>=^<2*'d<V'*>> (13) 

and 

TV- 7 



j*0 

20 



^) = IC/tf(T-0. (14) 



25 



C. Model of the Receiver Front End 

[0029] Taking into account (3), (4) and (1 ), the received signal (2) can be expressed as 



E 

30 **0 

where 

35 



X$ * Z ^('M'-O + "to = Z bj* gt,T=t-iT) + *<f) (15) 



/V-7 



j(/,r) = c(/,r)®/<r) = J c„(f)/<r-rj (16) 



40 



represents the time-variant equivalent channel impulse response at the receiver antenna. Apart trom the random phase 
shift introduced by the channel itself, imperfection and instability of the transmitter and receiver oscillators may give 
rise to a carrier frequency offset v. Denoting by f(t) the equivalent impulse response of the receive filter, te the result 
of all the filtering stages in the receiver front-end, the received signal after down conversion to baseband can be then 
45 expressed as 



so 



j{t) = e**"">M9 t\t) = bij<4t,t = t-il) + lit) , (17) 

i 

where 6 is a constant phase offset, 
55 q(t,?)= f • p(f,x)<2> f(x)= (18) 
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n{t)=/ 2 " vT * e) .w(t)<S> f{t)= (19) 
and ei 2jivt ji q(t,t) is the equivalent impulse response at the output of the receiver front-end. 

5 

CHANNEL ESTIMATION 

A. Channel Estimation and Tracking for a GSM receiver 

10 [0030] If the frequency response of the receiver filter f(t) is approximately flat within the bandwidth of x(0, F(f) = F[f 
(t)] = /for l/l < IB/2+vt, from (17)-(19) we can write 



(20) 



where 

20 



25 



30 



far) = J ^(O^-O, r„(i)= c.W^. (21) 



[0031] In a GSM receiver, channel estimation and equalisation are conventionally performed on a derotated and 
sampled version of the signal r(t). Assuming Baud-rate sampling, and defining r(k) = f k r{kl) t n(k) = f k n{k7), one has 



35 

Letting 



40 



we obtain 

45 



50 



HM-fiLt-kT.T-tD (22) 



iW = e'*"" T £ 4 4*>*-') + (23) 



55 



[0032] The channel impulse response is usually estimated by means of correlative channel sounding, which employs 
the training sequence midamble of the GSM burst. If b n , n = 0, N Q -1 indicate the known training bits, the estimate 
of the €-th channel tap h(k,€) is given by 
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m = -jj- S b » -4* + ") - ^'" , " r A^.O + r(0 (24) 



where n 0 = ^ + (A/ 0 -1)/2, and 

70 

rW = ~^*„ -^+^). (25) 

is [0033] Once the channel response is estimated by training at the burst midamble, the variation of the channel taps 
within the burst is usually tracked by means of decision-directed algorithms, where tentative decisions at the equalizer 
output are fed back to update the initial channel estimate. As a result of the tracking procedure, the time variations of 
the channel are represented by the estimated equivalent response 

h(k,i) ~ J nvkT h(k,€) + y( *,<?)■ (26) 

ESTIMATION OF THE DOPPLER SHIFT 

25 A. Effect of the Doppler Shift on the Phase Variation of the Channel Taps 

[0034] Assume the knowledge of the overall channel impulse response at the time instants t c =/c 0 Tand t, = k-,T. 
Denoting by h(k Q ,€) and h(k 1 ,€) the corresponding two values of the €-th channel tape, we consider the quantity 

30 

Yikesk,,*) = JiJc^Jfi^e). (27) 
[0035] The phase of (27) is the phase variation of h(k t ) in the time t r t Q . Using (22) and (21 ) 

Week,,*) = / ' ~g(k Jy t) J' gik^t) = g(*,,t)g (ko,i) , (28) 

40 

with g(k€) = g(t = kT t x=€T). Then substituting (16) into (28) 

n^k,,*) = ZI cAWc/i^rijitr-TjpieT-T,). (29) 

In the particular case where one can approximately assume the presence of one dominant line-of-sight component c n 
(r) with delay x n , from (6) we may write 

k i J) -v n * 0 7 

r<W>^> ( ( V^>| 2 

ss arg{r(/c 0 ,/c 7 ,f)} = 2n f Dn (k, - k 0 )T. 

However, when the tap-gains c n (t) represent diffuse scatting (see (5)), the relationship between the phase of the RHS 
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of (29) and the Doppler frequency is difficult to analyse. In order to identity a suitable measure of the Doppler, we 
observe that, taking into account (11 ), 



n~0 



N-1 



= RXkJ-KT) £ C n >\f(eT-T n )\ 3 = R^T-KT) ©^7), (30) 



where 

75 



20 



R c ( t 1- t 0)= J 0( 2 * 'o< f 7-'0» (31) 



= <M0 0 | /(r) |' = £ | /<r - r J f . (32) 

1^0 



25 Therefore, we conclude that with the classical Doppler spectrum 

arg{ 3 HKr^t) } = arg{ i* tfft,*) ^ ( V) > = <? (33) 

It is worth noting that the validity of (33) depends on the symmetry of the Doppler spectrum, and not on its particular 
shape. 

[0036] The above discussion suggests that an estimate of the Doppler should be based on a measure of the statistics 
of the phase variations of (27), that is, of the frequency variations of the €-th channel tap. We define the Doppler 
55 estimator 

2;c# D (/c 7 - k 0 ) T = ( E{(arg{Q(r)}) 2 }) 1/2 , (34) 

40 where Q(T) indicates a generic function of H/r^k-, ,*),€= 0, L-1. A simple way to decrease the influence of the taps 
with lower amplitudes is to choose Q(T) as the sum of the quantities r (k&kpi) for all channel taps. This gives 



45 *=* 



-(*<«i(£rf4.4/<*.0> )*>)*'• < 35 > 



so 

i-0 



[0037] The relationship between y 0 in (35) and the maximum Doppler frequency f D is illustrated in Figure 2 for the 
55 GSM type urban (TU) channel model. The data have been collected from the results of several simulation runs. 

[0038] The consistency of the Doppler estimates (35) for different multipath delay profiles may be improved by taking 
the function Q(T) = T(k Qt k v l) with I = argmax { \(g{kj) I ). 
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[0039] The corresponding estimator results 



= ( ^( «b{ ^,1) ^(^,7) } )*> (36) 

[0040] The dependence of the estimated Doppler shift $ D on the maximum Doppler frequency f D for the TU channel 
is shown in Figure 3. 

B. Measure of the Doppler shift using the Estimated Channel Response 

[0041] A practical Doppler estimator based on the rule (35) should make use of the estimated channel impulse 
response. Approximating the statistical expectation with a time average by sliding window, the estimator (35) becomes 

20 

Zrvofk, -WTm ( J-X (arg{ £ Tj^t)) )* )** = 

= ( —if (arg{f An^.^ V^O} )* )**. (37) 



70 



IS 



30 where h^k&f.) and h m (k v f.) are the estimated channel responses at the edges of the m-th burst, and the quantity 

3S 

is averaged over /bursts. Now, from (26)-(28), in the presence of a carrier frequency error we have 

k^tyki^t) * e^"<*' ™ T g(*o,*) + a**,,*,,*) . (38) 

with 

In the case of high signal-to-ndise ratio, neglecting the noise term ^(/co./t,,^), (38) gives 



arg{|; k^iykiKJ)} * 2tu(A>-^)7*+ arg{£ (*„,£) } - (39) 



55 



Therefore the Doppler estimate (37) is biased by the presence of a carrier frequency offset v. To overcome this problem 
we would need an estimate of the carrier offset. In fact, if such as estimate v m is available for the m-th burst, the Doppl 



er 
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estimator can be modified to 

2tf D fc -k 0 )T~ ( -^f; (arg{I k ]t l)}~ 2rrv m (k J -k^T) 2 V = 

= ( -if: hJLM-K'iM)-*tM-WY y 2 - 

10 M t-o 

75 

which is equivalent to (35), provided that the carrier offset estimate is not biased by the presence of a Doppler shift. A 
technique to derive an estimate of the carrier frequency offset to be used in the above Doppler estimator is illustrated 
20 in the next subsection. 

[0042] The above discussion applies to the estimator (35). It is straightforward to derive the corresponding result for 
equation (36) 

2n$ D (k, - J^JT m ( -Iff ( ars^fA, *„/)> - 2w m (k, - k^T )' )** (41) 

30 

C. Estimation of the Carrier Frequency Offset 

[0043] Assuming £(C(/fo*j.t)} * 0 t from (38) follows that 

Therefore, recalling that from (33) in the case of Rayleigh fading with symmetric Doppler spectrum 

40 

arg{ > = 0 . < 42 > 

we have the carrier frequency offset estimator 

so 

55 
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^v(*,-W= arg{ HXrjlc,,*,,*)) } 



- arg{ ^-^g^^,/)^^^) } = 



= Zw^-k^T + arg{ ^SMff(KJ) }, (43) 



20 and, provided that (42) holds, we obtain an estimate of the frequency offset that is independent of the Doppler shift 
Zzt^-k^T- arg{ kK^ ^i^i)} } « i3rK* 7 -4,)7\ (44) 



25 



30 



[0044] A practical estimator is again obtained by substituting the ensembled averages with time averages, according 
to 



3rv(*,-^)r* arg{ ±± £ r OT (^,*„£) } = 



40 



35 

= 4tL £ A-(*.0 • V<4.4 >, (45) 

A A 

where hj^k^i) and h m (k v t) are estimated channel responses at the beginning and at the end of the ro-th burst, and 
the quantity 

is averaged over Mbursts. 

[0045] An algorithm of the kind of (45) is usually implemented in the digital receiver of a GSM base station to correct 
so the frequency offset due to the transmitter and receiver oscillators. We observer that the quality of the estimation mainly 
depends on the effectiveness of the operation of noise suppression by time average, and on the tracking performance 
of the decision-directed algorithm that provides the channel estimate at the edges of the burst. 

SIMULATION RESULTS 

55 

[0046] The performance of the estimators discussed above in the section ESTIMATION OF THE DOPPLER SHIFT 
have been assessed by computer simulations for the GSM transmission channel The simulator includes GMSK mod- 
ulator, time-varying mulipath channel, additive white Gaussian noise, and receiver front-end. The equivalent noise 
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power spectral density at the antenna connector has been set to -116 dBm. At the receiver, after the equivalent base- 
band model of the RF and IF filters and after demodulation affected by a frequency error v, the signal is quantized and 
fed to a trellis equalizer. For each burst, an estimate of the channel is first computed by means of the GSM training 
sequence midamble. The estimated channel is then used to produce estimates of the transmitted bits, which in turn 
5 are fed back to track the variations of the channel impulse response. The following results have been obtained em- 
ploying an LMS tracking algorithm. 

A. Carrier Offset Estimator 

io [0047] In agreement with the analysis presented above, the simulation results for the GSM TU channel model show 
that the carrier offset estimator (45) is substantially independent of the Doppler shift. Figures 4 and 5 illustrate the 
performance with ideal frequency hopping, for a received signal level of-84 dBm. Similar results are obtained without 
frequency hopping. The behaviour of the estimator is summarised in Figure 6 and 7. 

75 B. Doppler Shift Estimator 

[0048] The performance of the Doppler estimator (40) has been assessed for the case of the GSM TU channel. In 
the simulation, the estimation of the carrier offset has been obtained from (45). As shown in Figures 8 to 11 , the Doppler 
estimator is not biased by the presence of a carrier frequency offset. From Figures 10 and 11 , one can derive that, in 
20 the case of frequency hopping, the estimated Doppler may allow to separate low mobility (3 Kmyh) and high mobility 
(100 Km/h) users at a received signal level of-104 dBm. The results collected from the simulation of 10,000 bursts at 
different signal-to-noise are summarised in Figures 12 and 13. 

[0049] Similar results are obtained with the Doppler estimator (41 ). The performance with TU channel and frequency 
hopping are shown in Figures 14 and 15. 

2S 

Description of a Preferred Embodiment 

[0050] Figure 16 shows a part of a base station having an equalizer 32 and a channel estimator and tracker 34, a 
carrier offset estimator 36, and a mobile user allocator 38 as is conventional. The equalizer 32 and channel estimator 
30 34 both receive incoming signals, and the estimator is connected to provide a channel estimation signal both to the 
equalizer and to the carrier offset estimator. 

[0051] There is also a Doppler shift estimator 40 connected to receive signals from the channel estimator 34 and 
also from the carrier offset estimator 36. The output ol the Doppler shift estimator is connected to the allocator 38. 
[0052] Figure 17 indicates the basic steps of signal processing, considering one burst and/or time slot. The m-th 

35 burst is indicated, and the direction of increasing k where k is the bit number. The burst is received through a channel, 
which equivalent discrete-time model has L-1 channel taps, the taps being spaced by a time T The real and imaginary 
parts of each tap evolve in time independently of each other and of the other taps. Within a burst, by calculations from 
the central bit ko outwards to either end of the burst, bit k 1t the speed of the mobile can be estimated. The objective, 
as explained above, is that the phase variations, derived from the frequency variations of -th channel tap are to be 

40 measured. 

[0053] The first step is that the estimator 34 estimates the -th tap of the channel impulse response for the time indexes 
ko, k v namely h^t) and h(k v €). The quantity r^k^) is then computed in step 52, (equation 27). 
[0054] The function Qfnko.k-, X) is calculated in step 54 (equation 34), and its argument derived in step 56 providing 
the first signal. 

45 [0055] Simultaneously, the carrier frequency offset estimate is derived in blocks 58, 60 and 62. In step 58 the sum 
of the quantities r(k 0 ,k-,/) is computed, for € = 0, .... L-1 ; its time average is performed in step 60, and the argument 
of the resulting signal 9 is taken in step 62, (equation 45) providing the second signal. The estimated carrier offset at 
the output of block 62 is subtracted in a subtract junction 42 and the squared value of the difference (64) is time 
averaged (66) and fed to a square root device 68 providing the estimated (RMS) Doppler shift. The estimated speed 

50 of the mobile can then be derived from the estimated Doppler shift, and this information can be used in allocating the 
mobile user to a microceli or a macrocell as appropriate. 

[0056] It will be appreciated that while it is known how to determine an estimate of the carrier Irequency offset in a 
mobile system, it is not previously been exploited, the fact that under the assumption of symmetric Doppler spectrum, 
the offset is independent of Doppler shift. The present invention depends on this realisation. 
55 [0057] Several further modifications of the invention will be readily apparent to those skilled in the art, in the light of 
the foregoing disclosure. The scope of the present invention should be interpreted from the following claims, as such 
claims are read in the light of the disclosure. 
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Claims 

1 . A method of determining the speed of a mobile user in a mobile telecommunications system comprises the steps of: 

5 (a) deriving a first signal related to a phase variation, in a time interval which is a burst or fraction of a burst, 

of one or more taps of an equivalent channel impulse response at the input of a channel equalizer in a base 
station receiver which receives the signal, 

(b) deriving a second signal related to the carrier frequency offset of the signal received in the base station 
receiver, 

10 (c) subtracting the second signal from the first signal, and 

(d) from the difference, deriving a third signal related to the Doppler frequency shift. 

2. A method according to Claim 1, in which the first signal is derived for each time interval by computing for each 
channel tap a quantity, the phase of which is the phase difference experienced by the tap in that time interval; 
computing the sum of said quantities; and deriving the argument of said sum. 



75 



20 



25 



A method according to Claim 1 , in which the first signal is derived for each time interval by computing for each tap 
a quantity, the phase of which is the phase difference experienced by that tap in that time interval; selecting the 
quantity corresponding to the tap having the highest amplitude; and deriving the argument of that quantity 

A method according to any preceding claim, in which the second signal is derived for each time interval by com- 
puting for each tap a quantity, the phase of which is the phase difference experienced by that tap in that time 
interval; computing the sum of the quantities; computing the average of the sum over time intervals related to the 
last burst and several previous bursts; and deriving the argument of the sum. 

A method according to Claim 4, in which the third signal is derived from the difference between the first and second 
signals by squaring said difference; deriving the average of the squared value over time intervals related to the 
last burst and several previous bursts; and deriving the square root of the average. 

A method according to any preceding claim, in which the channel impulse response is estimated at both edges of 
each time interval. 

A method according to any preceding claim, in which a plurality of third signals are derived relating to the Doppler 
frequency shift of signals received by the base station receiver along different multipath channels. 

A mobile telecommunications system comprising at least one mobile transmitter/receiver and a base station having 
a transmitter/receiver, the base station receiver having a channel equalizer, there being further provided first means 
to derive a first signal related to a phase variation, in a time interval which is a burst or a fraction of a burst, of one 
or more taps of an equivalent channel impulse response at the input of the equalizer; 

second means to derive a second signal related to the carrier frequency offset of the signal received in the 
base station receiver, 

subtraction means to subtract the second signal from the first signal; 

and calculation means to calculate from the difference a third signal related to Doppler frequency shift of a 
45 signal received from a mobile transmitter/receiver. 

9. A system according to Claim 8, in which the calculation means is arranged to calculate the speed of the mobile user. 



30 6. 



35 



40 



SC 



10. A system according to Claim 8 or Claim 9, in which the calculation means is connected to allocation means in the 
base station by which a mobile user is aiiocated to a microcell or macrocell of the telecommunications system. 



55 
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